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ABSTRACT 

The dwarf stars in the 26 year period binary a Com were predicted to eclipse each other in early 
2015. That prediction was based on an orbit model made with over 600 astrometric observations using 
micrometers, speckle interferometry, and long baseline optical interferometry. Unfortunately, it has 
been realized recently that the position angle measurements for three of the observations from ~ 100 
years ago were in error by 180 degrees, which skewed the orbital fit. The eclipse was likely 2 months 
earlier than predicted, at which point the system was low on the horizon at sunrise. 


1. INTRODUCTION 

a Comae Berenices has long been suspected of ec lips- 
ing (|Struvelll875l iHartkoof etahll 19891 iHo LHei tlIT 09f it ) de¬ 
spite being a 26 year binary, due to the system having 
an inclination extremely clo se to edge-on . The most 
recent orbital calculation by iMuterspaugli et all (|2010T ) 
showed the eclipses to be highly likely, with a pre- 
d icted closest projecte d approach in late January 2015 
(jMuterspaugh & Henrvf2014h . As the event approached, 
three findings showed the prediction likely to be in er¬ 
ror. The source of the mistake was that among the over 
600 observations used to determine the orbital model 
by iMuterspaugh et all ((2010h . the measurements from 
1896.33, 1911.4, and 1937.16 were listed with position 
angles in error by 180 degrees. All 3 of these were the 
last measurements made before closest approach of the 
binary at their respective epochs, marking a transition 
from measurements with position angles near 192 de¬ 
grees to those with 12 degrees or vice-versa. The stars 
are nearly equal magnitude, making such mistakes pos¬ 
sible. While efforts were made to find and correct such 
errors by examining fit residuals (see FigureQ] top panel), 
these three were missed because the orbit model skewed 
to compensate, as is possible near closest approach, and 
thus the fit residuals escaped detection (several other 180 
degree discrepant measurements were successfully cor¬ 
rected through this method). The final orbital solution 
was similarly skewed, which caused errors in the timing. 

2. CATCHING THE MISTAKE 

Three findings led to identifying these erroneous mea¬ 
surements and the skewed orbit solution. First, Henrichs 
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Figure 1 . Fit residuals (observed-computed) for the positional 
angles of a Com for various orbit models and 180 degree ambigu¬ 
ities. The horizontal axis is modified julian day and the vertical 
axis is residual angle in degrees. In each case, the 3 problem¬ 
atic measurements from 1896.33, 1911.4, and 1937.16 are circled. 
( Top) Fit residuals as th e data were presented for the model by 
IMuterspaugh et al.l dug). Even though the three measurements 
were input with 180 degree errors, the model fitting procedure re¬ 
sulted in a model skewed enough that the measurements do not ap¬ 
pear as obvious outliers. (Middle) Fit residuals for a model fit using 
only the separation measurements and PHASES observations, then 
applied to the positional angle data. The problem observations are 
clear. (Bottom) After fixing the 180 degree ambiguities of the three 
measurements, residuals to the currently presented orbital model. 
Measurements from sources other than the PH ASES program and 
CHARA array are from the WDS IlMason et al.lf2001fi . 

and Wijngaarden (preprint, to be submitted) refit the 
observations obtained from the WDS ignoring all posi¬ 
tion angle measurements (except to choose a positive or 
negative sign for the separations) and refitting an orbital 
model to the separation-only measurements (with incli¬ 
nation set to exactly 90 degrees and ignoring the longi¬ 
tude of the ascending node). Ignoring the position angle 
measurements turned out to be a fortunate choice, since 
the three erroneous position angles no longer skewed the 
fit. In a private communication to Muterspaugh on New 
Year’s Eve 2014, they speculated that the eclipse had in 
fact already passed based on this separation-only orbit 
evaluation. However, their orbital solution did not agree 
well with a full fit including the position angles; seeding 
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a full fit with their parameters led to fitting iterations 
which eventu ally con ver ged o n the orbital solutions of 
the iMutersuaugh et all (12010T) model. For two weeks it 
was unclear how to resolve the discrepancy. 

Second, in early January 2015 the system was observed 
by the Navy Precision Optical Interferometer (NPOI, 
lArmstrong et al.l 1199 81. The visibility trends indicated 
a b inar y separati on much larger than anticipated by the 
iMuterspaugh et ahl (|2010t ) orbit model, and it was un¬ 
clear if the system was still in approach or had already 
passed conjunction. 

T hird, a separati on m easu rements by the CHARA ar¬ 
ray (Iten Brummelaar et al.l 120051 1 on 16 January 2015 
found a separation of 45.53 milli-arcseconds (mas) and 
position angle of 192.85 degrees (which could be 180 
degrees ambiguous, though this value is now consistent 
with our revised orbit below), at a time when the pre¬ 
dicted separation from the 2010 orbit was only 7 mas. 
After this manuscript was first drafted, two additional 
measurements were obtained by CHARA and have been 
included in the orbital solution. These also show the 
binary to be growing further apart. In the CHARA ob¬ 
servations, the fringe packet of component “B” is about 
5% larger than what has been historically defined by the 
vast majority of observers as “A”; this is a function of 
the stars’ differential brightness, diameters, and also may 
be impacted by the fact that the CHARA measurements 
were at infrared wavelengths. To avoid confusion, the 
historical designations are maintained. 

As a result of t hese finding s, the original data set used 
in IMuterspaugh et all (1201011 was refit with the uncer¬ 
tainties for position angle artificially increased to ab¬ 
surdly large values (to 10,000 degrees), essentially remov¬ 
ing them from the fit, and a model was evaluated using 
just the separation measurements (with no positive or 
negative signs) and the PHASES measurements (which 
were made at a time when the quadrant was unambigu¬ 
ous). The resulting model was used to calculate pre¬ 
dicted position angles fo_r_ all observati on ti mes for com¬ 
parison with those in IMuterspaugh et all ( 2010 b It was 
discovered that the measurements from 1896.33, 1911.4, 
and 1937.16 were listed with position angles in error by 
180 degrees, namely 21.1, 196.4, and 208.3 degrees. The 
corrected values are 201.1, 16.4, and 28.3, respectively 
(see Figure [l] second panel). A new fit was performed 
with the position angle uncertainties returned to their 
original values, resulting in an improved \ 2 goodness-of- 
fit. 


3. UPDATED ORBIT 

The new orbital parameters for a Com were obtained 
to a combined fit based on the (corrected) non-PHASES, 
PHASES, and new CHARA array measurements (MJD, 
p (arcseconds), 0 (degrees): 57038.4870368, 0.04553 ± 
0.00081, 192.85± 1.02; 57043.4431444, 0.04855T0.00065, 
193.78 ± 0.77; 57044.48408989, 0.04947 ± 0.00043, 

193.19 ± 0.50). The results are presented in Table [T] 
There are 1217 degrees of freedom and y 2 = 1110. For 
comparison with the previous model’s epoch of perias- 
tron passage, which was listed one full orbit prior, the 
new model predicts MJD 47615.3±3.1 (the increased un¬ 
certainty compared to Table [1] reflects the uncertainty of 
the period which impacts this). The leading indications 
of the eclipse timing being incorrect are the resulting de- 



Figure 2. (Left) Histogram of the time of the eclipse mid-point, 
versus days since 20 November 2014 (MJD 56981). (Middle) His¬ 
togram of the closest projected separation of the binary (maximum 
eclipse) in units of milli-arcseconds. (Right) Histogram of the du¬ 
ration of the eclipse from ingress at 0.7 mas separation to egress 
at the same, in units of hours. 

creases to both the period and the epoch of periastron 
passage. 

4. NEW ECLIPSE TIMINGS FOR 2014 AND BEYOND 

Based o n the u pdated orbi t, the method outlined in 
IMuterspaugh &; Henrvl (I2014D was repeated to predict 
when the system did or will eclipse. We created 100,000 
random sample sets of binary orbit parameters based 
on the orbital model in Table |T| In each set, random 
values were selected for each orbital element using a 
Gaussian-distributed random number generator with lcr 
width, corresponding to the parameter’s formal uncer¬ 
tainty and centered on the best-fit value (e.g., values of 
the period were selected as 9443.1 + 3.0 x g days, where g 
is a standard normal deviate random number). The re¬ 
sulting set of parameters was then used to calculate the 
sky-projected separation of the binary every minute from 
MJD 56955 to 57003 (48 days), to ensure ail likely eclipse 
times were included. For each set we recorded the time 
of closest approach, the distance of closest approach, and 
the duration over which the binary separation was less 
than 0.7 mas (the approximate diameters of the stars). 
If any set failed to produce a minimum separation less 
than 0.7 mas, it was flagged as non-eclipsing, though in 
zero trials did this happen. 

The most likely time of eclipse was MJD 56981 (20 
November 2014, UT). Unfortunately, the observing sea¬ 
son for TSU’s photometric observations of a Com did not 
begin until five days after this predicted time of eclipse 
due to the star’s high air mass near the eastern horizon 
at sunrise. The observations were acquired with TSU’s 
T4 0.75m Automatic Photoelectric Telescope (APT). T4 
successfully observed a Com in good conditions for six 
consecutive nights beginning 25 November 2014, UT. 
Those six observations scatter about their mean with a 
standard deviation of 0.0031 mag and show no evidence 
for dimming. APT observations beginning a week or so 
earlier would have been difficult but perhaps not impos¬ 
sible. 

The secondary eclipse is now also a possibility, and is 
predicted to occur just eleven short years to the week of 
when the error in the 2014/2015 timing was discovered. 
The 100,000 simulations were repeated for a 48 day win- 
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Table 1 

Visual Orbit Parameters 


Period (d) 

To (HMJD) 

Semimajor Axis (arcsec) 

Eccentricity 

Inclination (deg) 

uj (deg) 

D (deg) 

crp 

°T 0 

CTa 

CTe 


er oj 


9443.1 

57058.37 

0.67160 

0.51104 

90.0576 

100.578 

12.215 

(3.0) 

(0.51) 

(0.00038) 

(0.00069) 

(0.0094) 

(0.027) 

(0.014) 


Note. — The revised model parameters and fit uncertainties for the binary orbit of a Com. 



Figure 3. Same as as Figure \2 \ except for the secondary eclipse 
centered at MJD 61051. 




Figure 4. Same as as Figure [2] except for the eclipse centered at 
MJD 66424. 

dow beginning on MJD 61032. The secondary eclipse 
has only a 3.3% probability of happening- in only 3318 
simulations did the stars come within 0.7 mas of each 
other. If the secondary eclipse does occur, the mean pre¬ 
dicted time of eclipse is MJD 61051 (11 January 2026); 
see Figure [3] 

The next primary eclipse is in late September 2040, 
a time of year which makes this system quite difficult 


to observe from Earth. Cameras on distant spacecraft 
could be used instead. The 100,000 simulations were 
repeated for a 48 day window beginning on MJD 66398. 
The average time of eclipse is MJD 66424 (27 September 
2040); see figure 01 

Figures [2j[3l and [4] show histograms of the results from 
100,000 simulated orbital parameter sets for all three 
events in terms of the time of eclipse minimum, the min¬ 
imum projected star separation (maximum eclipse), and 
the duration of the eclipse. 

We would like to thank and apologize to the many am¬ 
ateur and professional astronomers who were monitor¬ 
ing the system in anticipation of the eclipse. We thank 
Chris Farrington, Gail Schaefer, and the CHARA team 
for providing a timely measurement of the binary separa¬ 
tion. We thank Bob Zavala, Christian Hummel, and the 
NPOI team for observations helping to lead to these find¬ 
ings. We thank Brian Mason, Henrique Schmitt, Simon 
Albrecht, and Francis C. Fekel for reading early versions 
of this manuscript. 
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